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The -B* mass and lifetime measurements using the decay B^r — > J/ipn at the LHCb 
experiment were studied. About 310 signal events are expected for a data set which 
Ch ■ corresponds to an integrated luminosity of 1 fb _1 , with a B/S ratio equal to 2. Based 

on these data, the B^r mass and lifetime can be measured with expected statistical 



errors below 2 MeV/c and 30 fs, respectively. 

1 Introduction 



The meson is the ground state of the meson family formed by two different heavy flavor 
I ■ quarks, the anti-6 quark and the c quark. The mass of the Bj meson has been estimated 

Q-ij theoretically to be in a range of 6.2-6.4 GeV/c 2 [2J. Recently unqucnched lattice QCD [3] 

gave the most accurate prediction of M(B+) = 6278(6)(4) MeV/c 2 , where the first error is 



statistical, and the second is systematic. The be mesons are open-flavored and the excited 
states below the BD threshold can decay only through the electromagnetic or hadronic 
interactions into the ground state, Bj, which can decay only weakly and has a relatively 
f--*) ■ long lifetime. The estimated B+ lifetime is in a range of 0.4-0.8 ps [5J. At present, the best 

£ — \ prediction is given by QCD sum rules [3J: t(B+) = 0.48 ± 0.05 ps. Since the be mesons 

CO ■ carry two different heavy flavor quarks, their production [5] is more difficult than that of 

the heavy quarkonia. Including the contributions from the excited states, the cross section 
of the B+ meson at LHC was estimated [5] to be at a level of 1 fib, and is one order higher 
than that at Tevatron. This means that O(10 9 ) mesons can be anticipated with 1 fb _1 
of data at LHC, which is sufficient to study the be meson family systematically. 

The B+ meson was observed in the semilcptonic decay modes i? + — > J/ip([i + fi~)£ + X {I = 
e, /i) by CDF [B] at Tevatron. The measured mass and lifetime are consistent with the 
theoretical predictions. More precise measurements of the B+ mass have been performed 
recently using the fully reconstructed decay B+ — > J/ijj(fi + n~)ir + by CDF [7] and DO jS], 
giving M(B+) = 6275.6 ± 2.9(stat.) ± 2.5(syst.) MeV/c 2 and M(B+) = 6300 ± 14(stat.) ± 
5(syst.)MeV/c 2 respectively. The B+ lifetime measurements using the semilcptonic decay 
were also updated with more data, giving 0.448^0 036( s t a t-) ± 0.032(sys.) ps by DO [9] and 
0.475j;g;gg(stat.) ± 0.018(sys.) ps by CDF [Hj]. 

In this paper, we report our performance study of the B+ mass and lifetime measure- 
ments using the exclusive reconstruction of the decay £?+ — > J/iJj(h + fi~)ir + at LHCb. The 
LHCb detector is described elsewhere [11], hereafter is a summary of the LHCb detector 
performances which are relevant to this study. Based on Monte Carlo studies, the relative 
momentum resolution Sp/p was found to be 0.35%-0.55%, depending on the momentum. The 
primary vertex can be measured with a precision of about 10 /im in the plane transverse to 
the beam line, and of about 60 /im in the beam direction. The average resolution of the im- 
pact parameter respect to the primary vertex can be parameterized as <7ip = [14+ (35/pr)] 



1 Charge conjugate states are implied throughout this paper, if not specified differently. 
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fj,m where pt is the transverse momentum respect to the beam direction expressed in GeV/c. 
The muon identification efficiency was measured to be e(p — > fi) =94% with a misidentifi- 
cation probability e(ir — > fi) =3%. 

2 Event selection 

The simulation of the pp collision and the subsequent hadronization were performed with 
Pythia [IS]. £?+ meson was generated with a dedicated B c generator, BCVEGPY (T3] . 
The decay of any b hadron produced was done by EvtGen [14] . The response of the detector 
was simulated by Geant4 [15]. In this study, the specific b — > J/^p(p + fi~)X (including 
Bd — > J/4>X 7 B s -> J/ipX, B u — > J/ipX, Ab -> J/ipX), the inclusive J/-ip{n + H~) and inclu- 
sive bb MC data samples were simulated for the background study. 

The event selection [16] began with selection of good tracks based on the y 2 of the track 
fit, and requiring a pt greater than 0.5 GeV/c. A pair of two tracks with opposite charges 
and identified as p were used to reconstruct the J/ip meson. J/ip candidates were required 
to have an invariant mass in the range from 3.04 to 3.14 GeV/c 2 and a vertex fit quality of 
Y 2 /ndf < 9. Then a track identified as pion from the rest of the event was combined with 
the J/tp candidate to reconstruct the B+ meson with a vertex fit cut of x 2 / n df < 12. The 
impact parameter significance of the reconstructed B+ track with respect to the production 
vertex must be less than 3. For events with more than one primary vertex, the production 
vertex is given by the primary vertex to which the B^ candidate has the smallest impact 
parameter. To suppress combinatorial background further, pt(pl) was required to be greater 
than 1.0 GeV/c, pt(tt) and pt(B+) were required to be greater than 1.6 GeV/c and 5.0 
GeV/c respectively, and the impact parameter significance of the 7r and J/ip mesons were 
required to be greater than 3.0 and 3.5 respectively. 

The total reconstruction efficiency for the signal, including the trigger efficiency, was 
found to bc (1.01 ± 0.02)%. In the M(B+)±80 McV/c 2 (about ±3cr) mass window, which 
was regarded as the signal region, the background to signal ratio B/S with 90% confidence 
level is 1 to 2. The upper limit of B/S was used in the following study. One can expect 
about 310 signal events with 1 fb _1 of data assuming the B+ total production cross section 
a(B+) = 0.4 fib and BR(S+ J/?/>7r+)= 1.3 x 10~ 3 . 

3 Mass measurement 

To improve the B+ mass resolution, the J/ip mass constraint vertex fit was used in the 
mass measurement. M(B ( t)±300 MeV/c 2 mass window was used in the following study. 
The un-binncd maximum likelihood method was implemented to extract the B+ mass. The 
invariant mass distribution of the signal was modeled by a single Gaussian, that of the 
background was parameterized as a first order polynomial. 

For the fit the signal events were taken from the signal data sample with the full Monte 
Carlo simulation and the background events were generated by a standalone simulation using 
the invariant mass distribution obtained from the full Monte Carlo simulation. Both the 
signal and the background were normalized to 1 fb" 1 of data. The fit gives the B+ mass to 
be M(B^) = 6399.6 ± 1.7 (stat.) MeV/c 2 , which is in good agreement with the input value 
of 6400 MeV/c 2 . The result of the fit is shown in Figured] 
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Figure 1: The candidate mass distribution together with the result of the fit described 
in the text. 

4 Lifetime measurement 

The J/ip mass constraint vertex fit biases the B^ decay vertex position because of the final 
state radiation in the J/ip decay and the energy loss of particles. As a consequence, the J/ip 
mass constraint vertex fit would bias the lifetime and was not used in the lifetime study. 

A combined mass-lifetime fit was used to measure the B+ lifetime. The invariant mass 
distributions for the signal and the background were obtained as in the mass measurement. 
The proper time (t) distribution of the signal was described by an exponential function 
smeared by resolution and multiplied by an acceptance function e{t) describing the distortion 
of t distribution caused by the trigger and offline event selection. The function form was 
derived from an independent Monte Carlo sample. The t distribution of the background was 
parameterized by the "KEYS" (Kernel Estimating Your Shapes) method [17]. The proper 
time resolution {at) distributions of the signal and the background were also considered to 
avoid the potential bias caused by the so called Punzi effects [18] , The fitting framework 
was tested with a standalone simulation program. The pull of the _B+ lifetime was found 
to follow a Gaussian distribution with fi = 0.004 ± 0.010 and a = 1.0122 ± 0.0072, which is 
consistent with the Normal distribution. 

As shown in Section [TJ the resolution of the impact parameter depends on the transverse 
momentum. Therefore, the impact parameter significance and the proper time acceptance 
e{t) depend on the pt{Bj) distribution. In order to evaluate the systematics introduced 
by this, a fit to which was generated with the pt distribution predicted for B + by 
Pythia, was made using the acceptance function obtained from the B+ data generated 
by BCVEGPY. Note that the Pt{B+) spectrum predicted by BCVEGPY is harder than 
Pt{B + ) spectrum predicted by Pythia. A standalone simulation study shows that in this 
case the B+ lifetime would be biased by about 0.023 ps. 

To reduce the dependence of the lifetime measurement on the p^{B+) distribution, the 
pr{B+) was divided into two intervals: 5-12 GcV/c and > 12 GeV/c. The event selection 
was re-optimized and for the data set with pt{B+) > 12 GeV/c, the impact parameter 
significance cuts are changed to > 2.0, > 2.5 and < 4.0 for pions, J/ip and £?+ respectively. 
The event selection results of the two pt intervals are summarized in Table [TJ The combined 
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Pt intervals of £> + 


5-12 GeV/c 


> 12 GcV/c 


Total efficiency £ tot (including trigger) 


(0.34 ±0.01) % 


(0.86 ± 0.02) % 


Signal yield (1 fb" 1 @ 14 TeV) 


100 


260 


B/S @ 90% CL 


3 to 6 


0.6 to 1.2 



Table 1: Summary of the event selection results in the two pr(B+) intervals. The total 
efficiency £tot in the two pt intervals are both defined with respect to the total number of 
the events in the full px range. 



mass-lifetime fit applied simultaneously for the two pt(B^) intervals reduces the bias in the 
lifetime measurement to about 0.004 ps. 

The fit to the 1 ftr 1 of data gives B+ lifetime to be r(B+) = 0.438 ± 0.027 (stat.) ps, 
which is in good agreement with the input value of 0.46 ps. The results of the fit in the two 
Pt(B^) intervals are shown in Figure [2J 





Figure 2: The results of the fit in the two pt(B+) intervals. The left column shows the 
invariant mass distributions, the right column shows the proper time distributions in the 
signal region. 
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5 Summary 



In summary, the event selection of the decay — > J/ip(/i + h~)it + at LHCb was studied, 
about 310 signal events are expected from 1 fb _1 of data, with a B/S of 2. Based on the 
selected data sample, the B~£ mass can be measured with an expected statistical error to be 
below 2 MeV/c 2 . For the B+ lifetime measurement, to reduce a potential bias caused by the 
limited knowledge of the pr(B+) distribution, the data sample was divided into two Pt(B+) 
intervals. The combined mass-lifetime fit performed in these two pt intervals simultaneously 
gives a statistical uncertainty of the lifetime to be below 30 fs. 
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